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oo : We investigate the depinning of a massive elastic manifold with d internal dimensions, embedded 

in a (d + n)-dimensional space, and subject to an isotropic pinning potential V(u) = V(|u|). The 
tunneling process is driven by a small external force F. We find the zero temperature and high 
temperature instantons and show that for the case 1 < d < 6 the problem exhibits a sharp transition 
Oh' from quantum to classical behavior: At low temperatures T < T c the Euclidean action is constant 

up to exponentially small corrections, while for T > T c , SEuci(d,T)/h — U(d)/T. The results 
are universal and do not depend on the detailed shape of the trapping potential V(u). Possible 
' applications of the problem to the depinning of vortices in high-T c superconductors and nucleation 

in d-dimensional phase transitions are discussed. In addition, we determine the high-temperature 
asymptotics of the preexponential factor for the (1 + l)-dimensional problem. 
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I. INTRODUCTION 

i-^j \ The decay of metastable states is a basic phenomenon of great generalitvtl with numerous applications in a large 
number of contexts, ranging from the -decay of false vacua in field theoryB, e.g. in cosmologyli, to the creep type 
motion of topological defects in solidarEl. At a given temperature T the inverse lifetime T of a metastable system can 
be written in the form T = Ae~ SEud ( T ^ h , with SeuciC? 1 ) the Euclidean action of the saddle-point configuration and A 
the prefactor determined by the associated fluctuations. In this paper we present a calculation of the full temperature 
dependence of the Euclidean action of a metastable d-dimensional clastic manifold pinned by an isotropic potential 
V(u) = V(|u|) in the presence of a small external force. Typical examples of physical realizations of this model 
system are a string pinned by a linear or columnar defect ((1+1)-, and (l+2)-dimensional problems respectively) or 
a membrane pinned by a plane ((2+l)-dimensional problem). 

With increasing temperature, the decay changes its nature from quantum to classical. The function Seuci(T) might 
be either a smooth function of temperature or exhibit a sharp kink with a discontinuity in its first derivative at the 
qq ' crossover temperature T c . The former case is called a "second-order" transition from quantum to classical behaviom, 
Q"\ , while in the latter case a first-order transition takes place. The wording "transition" is appropriate as the crossovers 
have all the features of mean-field phase transitions upon identifying the Euclidean action with the free energy. n 
The possibility of a first-order transition in a tunneling problem has first been discussed by Lifshitz and KaganQ 
in the context of quantum phase transitions in 4 He systems. The general theory and an appropriate criterion_for a 
first-order transition in ID Hamiltonian systems has been developed by MeshkovS and by Ioselevich and RashbaB, see 
also the paper of ChudnovskyEj. For the case of non-Hamiltonian systems this problem has been considered in Ref. 
|i~T| . Recently, the possibility of observing a first-order transition from quantum to classical behavior in spin systems 
has been discussed in Ref. [l^. ._. ._. 

Systems with many degrees of freedom have been studied by Garrigalij (see also Ref.EJ) and he finds that, for the 
case of tunneling of massive d — 2 and d = 3 elastic manifolds between two minima disbalanced by the action of a 
small external force, the Euclidean action-exhibits a discontinuous derivative, while for a d = 1 dimensional string 
the problem exhibits a smooth crossovenla. A natural counterpart of the bistability problem studied by Coleman, 
Garriga, and others is the generic problem of weak metastability, see Fig. 1. In this paper we then study the problem 
where an elastic manifold is pinned by an isotropic potential well and is driven by a weak external force. After the 
tunneling process (or thermal activation at high temperatures) the manifold becomes free, while in the previous case 
it is repinned again in the second valley. The approximation introduced by ColemanB and used by Garriga is usually 
termed the "thin wall approximation" , while in the present problem it is more appropriate to denote the technique 
used as a "thick wall approximation" . 

We find that at least for the case 1 < d < 6 the problem exhibits a sharp transition from quantum to classical 
behavior: The function SEuc\{d, T) is a constant (up to exponentially small corrections) at temperatures lower than 
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the critical temperature T c , while for T > T c , SEuci{dJT)/h — U(d)/T, with U the activation energy. The (1+1)- 
dimensional problem has been considered by SkvortsovEj. In the present paper we generalize the results of this work 
for an arbitrary number of internal degrees of freedom d and improve on the calculation of the preexponential factor 
for the (1 + l)-problem as the approximations made in Ref. [lj] turn out to be too rough to reach the correct result. 
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Fig.l The two generic decay problems of bistability (a) and weak metastability (b). 



The paper is organized as follows: in Sec. II A we formulate the problem. In Sec. II B and II C the zero temperature 
and thermal instantons are found and the corresponding Euclidean actions are calculated. In Sec. Ill the instability 



temperature of the thermal instanton is calculated and the full temperature dependence of the Euclidean action is 
found. Sec. IV is devoted to the calculation of the high-temperature asymptotics of the preexponential factor for the 
(1 + l)-problem. We summarize the results in Sec. [v] and discuss their possible applications to the thermal depinning 
of vortices in high-T c superconductors and nucleation phenomena in phase transitions. 



II. EUCLIDEAN ACTION AT LOW AND HIGH TEMPERATURES 

A. General expression 

Consider a massive manifold with d internal dimensions embedded in a (d + n)-dimensional space. The (real time) 
Lagrangian of the manifold can be written in the form 

n K J 

Here, u e TZ n is the (transverse) displacement field of the manifold, r € lZ d is the vector characterizing the internal 
degrees of freedom, p and e are the mass density and the elasticity of the manifold, respectively, and V(|u|) is an 0(n)- 
invariant trapping potential. The function V(u) is supposed to be monotonously increasing, V(0) — 0, V(oo) = Vq. 
The integration in Eq. (Q) goes over the volume fl of the manifold. Finally, F = (F, 0, . . . , 0) is the external driving 
force which we assume to be small, F <C Vq/cl, with a the characteristic radius of the pinning potential V(u). Specific 
realizations of this model system are strings moving in a plane ((1 + l)-problem) or in space ((1 + 2)-problem), or 
membranes pinned by planar boundaries or interfaces in 3D space ((2 + l)-problems). 

After the transformation t — ► —it and S — > — iS^uci the Euclidean action of the manifold reads 

+h/2T (22 •) 

s e M*,t)}= I */*< f(£) +V(M)-F.u • (2) 

-n/2T n ' 

After introducing the new variables x' = x/y^e and r' — r/^fp, the Euclidean action (|2|) can be rewritten in the form 
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S Eucl [u(x', r')] = p 1 / V/ 2 / dr> J dx' + \ + V(\n\) - P • u 1 . (3) 

Below we shall work with this Euclidean action and find its zero- and finite temperature instantons. 

B. Zero temperature instanton 

Let us calculate the Euclidean action corresponding to the extremal trajectories at a given temperature T. The 
variation of the action (j3|) yields the saddle-point equation 

o + A'u = F. (4 

We need to find the solution satisfying the boundary condition u(x', 0) = u(x', h/y/pT). The solution of Eq. (|J) can 
be written in the form u = (u, . . . , 0) . At zero temperature the instanton is spherically symmetric and introducing 
the new variable r 2 = t' 2 + r' 2 we obtain an ordinary differential equation for the function u(r), 

l(rV)' = u " + ^u' = ^-F. (5) 

The function u(r) must be non-zero and continuous, have a continuous derivative u'(r), and satisfy the boundary 
condition 

u(r), u'(r) — > 0, r — > oo. (6) 

The corresponding solution consists of two parts: a "macroscopic" solution for r ^ Rd (Rd is a parameter which will 
be determined later) and a "microscopic" solution for r > Rd- The former corresponds to displacements u > a (we 
remind the reader that a is the characteristic radius of the pinning potential), while the latter describes the bounce 
in the region u < a, where the potential is relevant. 

Consider the "microscopic" solution where we can neglect the (small) force F in the equation of motion (||) as 
compared to the potential term d u V . Furthermore, in the limit of small F the radius Rd of the nucleus is large (we 
will check this assumption at the end), such that we can neglect the term du' jr in Eq. (||), and the remaining equation 
takes the form 

Eq. (Q) is identical to Newton's second law. Accounting for the boundary condition ^) we obtain 

- V(u) = 0, (8) 

which is equivalent to the law of energy conservation for the ID conservative motion in classical mechanics. As 
V(u ^ a) ~ Vo, we can obtain the boundary condition for the function u(r) 

u'{R d ) - -V2Vo- (9) 

A simple integration then provides the microscopic part of the solution, once the parameter Rd is known. 

Next, we solve Eq. (||) in the region r < Rd with the boundary conditions given by Eq. (||) and u{Rd) ~ a. As 
we will see below, u(0) cx l/F, allowing us to use the condition u(Rd) — without changing the main term in the 
asymptotics of the Euclidean action in the limit F — > 0. In the region r < Rd we can neglect the term dV /du in 
Eq. (|5|) and the general solution can be written in the form 

u[r) = ~W^) +Cl + C2 - f{rl (10) 

where 
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./■(,■) = < H-l' d?61> (11) 

lnr, d = 1. 

The term C2/ (r) is always singular at the point r = 0, hence C2 = 0. Taking into account the boundary condition 
@ we obtain 

^.VH^+l) (12) 

and the quantum instanton can be written in the form 

»M = ^(^-0. (13) 

Substituting the function it(r) into Eq. ([|) and taking into account that r 2 = r' 2 +x' 2 we obtain the zero temperature 
Euclidean action 

5 E uci(T = 0, d) = — (d + 1) — — Vp^ot7 + O ttt , (14) 



(d + 3) \ F J y ^ F \F d 

where v4 d = 27r rf / 2 /r(d/2) is the surface area of a unit sphere in d-dimensional space. 

Neglecting the terms F or dV/du in Eq. (^) might seem questionable as there is an interval with the length of order 



Uf where dV/du 



F. However, if this interval is much smaller than Rd the approximation used above is valid: 



Assume that the potential V(u) at large u behaves as V(u) ~ Vq — B/u a , a > 0. The equation V' — F, F — > can 
be easily solved and yields the result up ~ (aB / F) 1 ^ a+1 , i.e., in the limit F — * 0, tif -C i?^ and the approximation 
made above is asymptotically correct. 



C. Thermal instanton 



Next, let us calculate the activation energy. The saddle-point solution is time-independent in this case, i.e., we 
can eliminate the term d 2 u/dr /2 in Eq. (|J). The remaining equation is identical to the (d — l)-dimensional zero 
temperature problem considered above and we can use the instanton (|l^) with d — > d — 1. We substitute this solution 
into Eq. (^) and obtain the expression for the activation energy, 

„ w= &-.f«)V D+0 av (i5, 



(d + 2) V F J \ Fd 

We remind the reader that for the thermal instanton Se uc \(T ^> T c , d) = HU(d)/T. 

Next, let us calculate the two characteristic temperatures T c (d) and T*(d). The former corresponds to the value 
of T, where the quantum and thermal exponents SEucl(d)/7i and U{d)/T match up, while the latter defines the limit 
of applicability of the zero temperature instanton (see Fig. 2) : at temperatures T > T* (d) the instanton Jl3| ) does 
not satisfy the periodic boundary condition u(x, 0) = u(x,h/T). Note that the thermal instanton is always a valid 
solution as it is time-independent and satisfies the periodic boundary condition identically. 

The exponents S-E UC \(d)/h and U(d)/T are equal at the temperature T c (d) given by the expression (we make use of 
Eqs. (H|) and @) 



27rd(d + 2) Vd+1 



Formally, the quantum instanton found in Section [I B is valid only at zero temperature, however, we can use it 



also at T ^ as long as the periodic boundary conditions are satisfied asymptotically. The zero temperature bounce 
solution then can be still applied if the periodicity in imaginary time r', h/ \fp T, is larger than the diameter of the 
nucleus 2R d , i.e., if R(d) < h/2 y /pT (see Fig. 2). 
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Fig. 2 Nucleus corresponding to the low temperature decay process. As long as the diameter 2Rd remains smaller than the periodicity 
in time h/^/pT (solid lines) the instanton (^) asymptotically satisfies the boundary conditions. At temperatures T > T*(d) = S/2 v /pi? ( i 
the boundary conditions are not fulfilled and the distortion of the nucleus' walls becomes relevant. 

Making use of Eq. ( |l2| ) we obtain the following expression for the temperature T* (d) at which one needs to take 
into account the distortion of the quantum instanton due to the periodic boundary conditions, 

T " w = 57^^' (17) 



III. INSTABILITY TEMPERATURE AND EUCLIDEAN ACTION 

In this section we calculate the instability temperature To(d) of the thermal instanton. At this temperature the 
saddle-point solution becomes time-dependent. Then we show that To(d) < T c (d) < T*(d) for 1 < d < 6 and, 
consequently, the problem exhibits a first-order transition from quantum to classical behavior. 

Let us expand the Euclidean action around the thermal saddle-point solution. As the bounce satisfies the Euler- 
Lagrange equations, the term proportional to Su (Su describes the distortion of the saddle-point configuration) is 
equal to zero. To second order in Su we obtain 

+h/2^pT 

SeuciIu^x') + Su] ~ 5 Eu cl[utfc(x')] + i J dr' J dx'Sm (ill - P^y^j Su i 

-h/2^pT Q' 

+h/2^pT 

+ ~ J dr' J dx'5u ± (h x - P^j Su ±1 (18) 

-h/2^/pT Q' 

where 8u\ — Su ■ (1, ... 0) and Su± = Su — 8u\ • (1 ... 0). The operator Hi can be written in the form 

d 2 V 



= -A' 



du 2 



(19) 



The instability temperature is determined by its only negative eigenvalue. The operator H± has only positive eigen- 
values. 

Let us calculate the lowest eigenvalue — cu 2 of the operator (^9|). Here, we follow the procedure described in detail 
in Refs. [l3 15, 17]. First, we show that the operator ( |l9| ) has d zero eigenvalues: the d functions <fri — duth/dx^, 



i = 1, ... ,d, satisfy the equation Hx4>i = 0. Due to the spherical symmetry of the instanton and, consequently, of 
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the "potential" in Eq. fll9|), the eigenfunctions of the d-dimensional Schrodinger operator (|19|) can be represented 
in the form = = dfUthd x '?, with ^(f) and li(f2) the radial and angular components, respectively, 

/2 



Ed 
i— J- 



For the d eigenfunctions corresp ondin g to the zero eigenvalue ^E'i(r) = ^(r) — dfUth °c f for f < R d -i 



(see Eq. ( |13| ) and the first paragraph in section II C ). In the region f > the eigenfunctions of the lowest levels 



are asymptotically (F — > 0) equal, allowing us to find the boundary condition for the lowest eigenvalue — to 2 from 
that of the zero eigenvalue wavefunction. The condition for binding the radial parts of the wavefunctions at the point 
r = Rd-i takes the form 



1 



f=R d - 



Rd-l 



The general solution of the equation 



in the region f < R d can be written in the form 



*-i(f) = f {2 - d)/2 [C 1 / (d _ 2)/2 (^f) + C 2 K {d _ 2)/2 {^f)} 



(20) 



(21) 



(22) 



where I u (z) and K v (z) are Bessel functions of imaginary argument. The function Ki d _ 2 y 2 (z) is singular at the point 
z = and hence C 2 = 0. (For the case d = 1 this analysis is not valid, however it can be easily shown that the equation 
for the lowest eigenvalue (see Eq. (|23|)) remains the same). Using Eq. pQ| ) we obtain a transcendental equation for 
the negative eigenvalue — uj 2 (d) = — \-n 2 pT^ (d)/h 2 , 



(d-2)/2 



and 



I(d- 2 )/ 2 {Hd) 
Md) 



-, fj, d = uj(d)R d -i, 



I'd 



hF 



2^nd \fpVo' 

Using Eqs. © and (0) and the definition of the parameter /i d , see Eq. (|23J), we find 

2 (d+l)(d + 3)/' d Vr(^±i) 



and 



/?r (d + 2)d 

Tg(d) _ d + 1 (jb d 
T c (d) d nxd 



(23) 



(24) 



(25) 



(26) 



TABLE I. Ratios T c (d)/T*(d) and T (d)/T c (d) for 1 < d < 6. For T (d) < T e (d) < T*(d) the system exhibits a sharp 
transition from quantum to classical behavior. For d = 7, both ratios exceed unity and no conclusion can be drawn regarding 
the order of the transition. 



dimensionality, d 


l-i d 


Xd = T c (d)/T'(d) 


r (d)/r c (d) 


1 


1.1997 


4- « 0.8488 

37T 


0.8998 


2 


1.6083 


| « 0.8333 


0.9215 


3 


1.9150 


& « 0.8594 


0.9457 


4 


2.1725 


iif « 0.8960 


0.9647 


5 


2.3993 


5°°° « 0.9357 

17017T 


0.9795 


6 


2.6048 


i - so "' 9759 

33614 U.t7iu» 


0.9918 



G 



In Table | we summarize the numerical values of the parameters Xd and To(d)/T c (d). We see that for 1 < 
d < 6 the relation To(d) < T c (d) < T*(d) holds and, consequently, the Euclidean action has the form plot- 
ted in Fig. 3: It is equal to 5 Eu ci(d) (see Eq. (Q) for T < T c (d) (see Eq. @), whereas for T > T c (d) 
we have <Seuci(^, T) — TiU{d)/T. The jump from the zero-temperature to the high-temperature instanton takes 
place at T — T c : the system exhibits a sharp transition from quantum to classical behavior, see Fig. 3. 




To(d) T c (d) 



T*(d) 



Fig.3 The Euclidean action as a function of temperature for 1 < d < 6. For T < T c the Euclidean action is constant and given by 
Eq- (0) ( see E q. (0) f° r an expression for T c ). For T > T c the action is equal to HU(d)/T, with U(d) given by Eq. (O). 



IV. PREEXPONENTIAL FACTOR FOR THE (l+l)-PROBLEM 



In this section we calculate the preexponential factor for the (1 + l)-problem in the-high temperature region, 
improving on the results of Ref. Il6|. We start from Langer's expression for the decay ratals (see also Ref. Il9]) 



TiT = 2T (1)- 



ImZ 



(27) 



where T (l) f=irE<$(d — 1), see Eq. ( p4| ) and Table || and Z is the partition function of the system. For a (l+l)-string 
we can writeli3EJ 



ar 



dx 



du 



Hi 



dx 



1/2 



det (S 2 S/5u 2 \ u=0 ) 



det' {5 2 S/5u 2 \ u=Uth{r) ) 



1/2 



exp 



T 



(28) 



where L is the length of the string and the prime indicates that we exclude the zero eigenvalue of the operator 
Hi — S 2 S/Su 2 \ u=Uth ( x y The activation energy U(l) is given by the expression 



U(l) = — V e ^o-p 



(29) 



and the eigenvalues of the operators H = 5 2 S/Su 2 \ u=0 = — dl + d 2 V\ u=0 and Hi = 5 2 S / 5u 2 \ u=Uth ( x ) = —d 2 + d 2 V\ Uth 
take the form 



Ao.r 



/2ttT \ 2 
A ,«+p(— n\ , n = 0,±l,±2. 



and 



/2ttT 



Ai, an = Ai >Q + py— — nj , n = 0,±1,±2..., 
respectively. Carrying out the product over n wc obtain the high temperature result 



(30) 



(31) 
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nr = 



T (1)L 
V2^T 



+00 



dx 



duth 
dx 



-1 1/2 






deti? 




det'Fi 



1/2 



exp 



U*(l) 



(32) 



where 



U*(l) = 17(1) 



dfc 6p(k) v k + efc 2 



(33) 



is the quantum renormalized activation energy, see Ref. |T^. In Eq. ( p3|) k = V"(0), <5/o(fc) is the difference in the 
continuous part of the spectral densities of the operators i7o and i7i , and the sum is taken over the positive discrete 
eigenvalues of the operator Hi (note that the operator H has no discrete eigenvalues) . At large k* the correction to 
U(l) diverges as Ink*, i.e., we have to-introduce a cutoff k* into the problettu-The determinant ratio in Eq. ([52]) has 
been calculated by Kramer and KulictLZl using the Gelfand-Yaglom formulaE!rE3 with the result 



dctiJo 



det'iJi 



F {V2kVo 
— exp 



(34) 



The coefficient of proportionality depends on the detailed form of the pinning potential V(u). The final answer for 
the decay rate then takes the form 



FL 




exp 



V^kVo U*(l) 



F 



T 



accounting for gaussian fluctuations we obtain an exponential enhancement of the decay rate. Note that U* 
see Eqs. («3f) and (29). The result (p5[) is valid for temperatures T below the kink activation energy Ej~ 



FU*(1) 



(35) 



l/F, 



(36) 



(here we assumed that the quantum correction to the activation energy U(l) is small and k ~ Vq/o 1 ). At higher 
temperatures we need to take into account the renormalization of the free energy due to thermal fluctuations: The 
free energy per unit length of an elastic string trapped in a poteptial well coincides with the ground state energy of a 
massive particle in the equivalent quantum mechanical problemca, where the temperature plays the role of Planck's 
constant. If the temperature is small, the ground state energy is close to zero. At high temperatures, satisfying the 
condition Vg < T 2 /ea 2 (which is the condition for a weak potential well, see Ref. ^4|), the ground state energy increases 
substantially and becomes comparable with the potential well depth. In this case the thermal instanton found in Sec. || 
is not applicable and the problem is more conveniently solved in the quantum mechanical formulation: Since in the 
ID quantum particle trapped in an attractive potential always exhibits a bound state, the (1 + l)-dimensional string 
problem at F = does not undergo a depinning transition at any temperature. 



It is instructive to compare the resulLi 
between two nearly degenerate minimatS 



..to that obtained for the case of thermally activated motion of a string 
3, see Fig. 1(a), where the prefactor shows a power-like dependence on 
the external force F. The different dependences of the prefactors follow fronipthc different low energy spectra: In the 
present case we have ~ Vo/Fa low-lying eigenvalues of order F 2 /Vp (see Ref.cZl), the product P(F) of, which behaves 
as In P(F) ~ —1/F, while for the problem studied in Refs. [l5| , p5| , p6| only one small eigenvalue existstil. 

Finally, let us compare the result ( |35| ) with that obtained previously by SkvortsovcJ where the semiclassical ap- 
proximation for the eigenvalues of the operators Hq and Hi has been used. The derivation in Ref. [l(] differs from the 
present one in two respects: i) The calculation has been based on the low temperature expression 



HT = 2T- 



ImZ 



(37) 



for the decay rate, such that the final high temperature result lacks the correct classical limit (r oc 1/h, see Eq. (57) in 
Ref. |l6| ). ii) Correcting for the factor Tq/T, the remaining difference can be traced back to the use of the quasiclassical 
approximation in the calculation of the ratio (|35j), which produces a prefactor cx F 3 / 2 rather than the correct result 
cx F. Note that at zero temperature Eq. (|37]) is applicable, however, the semiclassical approach is not accurate enough 
to produce the correct expression for the decay rate (though it is still possible to obtain the correct exponential 
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enhancement of the preexponential factor and the correct quantum renormalized Euclidean action S*, see Ref. |l§ 
similar to the high temperature case discussed above). In principle, the complete prefactor in the low temperature 
quantum regime can be caclulated using the procedure suggested in Refs. |2^,|2{| The calculation of the determinant 
ratios can be reduced to the calculation of an infinite product of determinants of ID Schrodinger operators, which 
then can be calculated numerically. 



V. CONCLUSION 



As possible applications of the present problem we discuss i) the thermal depinning of vortices in high-T c super- 
conductors and ii) the nucleation of phase transitions in the vicinity of a boundary or an interface. 

i) The dynamics of vortices in high-temperature superconductors is dominated either by the dissipative or the Hall 
term. The vortex mass can be neglected. However at high temperatures the activation energy is independent of the 
dynamics (see Eq. ([l5])) and we can use Eq. ( |l5| ) for the description of the thermal depinning of a single vortex line 
(rf = 11 from a columnar defect (n — 2), or of a vortex sheet (d = 2) from an interface between superconducting 
graina23 (n — 1). We wish to point out that Eq. ([If]) can also be used for the activation energy of a single vortex 
depinning from a twin boundary. Note that a vortex pinned by a columnar defect ((1 + 2)-dimcnsional problem) 
again does not undergo a depinning transition at any finite temperature: the corresponding 2D quantum mechanical 
problem always has a bound state and hence the vortex remains localized at all temperatures provided that F = 0, 
see Ref. |2^. As our result is restricted to temperatures below the kink energy, our analysis is not in contradiction 
with this result. n 

ii) Consider the nucleation of a new phasdJ in the vicinity of a boundarytjor interface. Let us assume that there exist 
two competing (Z?-dimensional) phases, the stability of which are depending on the given external conditions. Initially, 
phase 1 is prepared, after which the external conditions are changed adiabatically, rendering phase 1 metastable. The 
nucleation of phase 2 takes a finite time and might be of two types: usual bulk nucleation or nucleation in the vicinity 
of the boundary. The latter can be described as the problem of depinning of a ((D — 1) + l)-dimensional manifold. 
The boundary between the two phases plays the role of the elastic manifold, while the difference in the chemical 
potentials plays the role of the small external force. It is interesting to study the competition between these two 
possible types of nucleation. The free energy Ud of the D-dimensional spherical nucleus of radius R is given by the 
expression 

U D {R) = eAoR - 1 - FV D R D , (38) 

where A D = 2tt d ^/T(D/2) and Vd = 2ir D / 2 /DT(D/2) are the surface and the volume of a unit sphere in D- 
dimensional space. Minimizing this energy with respect to R we obtain 



and 



2n D / 2 (D-lf- 1 (e\D-i 
UD = TiDj2) D <f) ' ^ 



This expression should be compared to Ud = U(D — 1) as given by Eq. ([Tq). After some algebra we obtain 

Ud = (D 2 -l)T((D-l)/2) f {D+1)/2 

U D 2(^+!)/ 2 D T{D/2) \V Q 



(41) 



It is reasonable to assume that e « Vq, i.e., the boundary tension e roughly matches up with the pinning potential 
Vq. In this case Ud/Ud > 1 for D = 2 and 3, i.e., nucleation at the boundary is more favorable than that in the 
bulk (the same result applies also for the weak pinning case Vq «C e). The analysis for the quantum case is more 
complicated: The kinetic energy of the nucleus can be written in the form f(R)R 2 , with f(R) a function of the 
nucleus' radius R. The process of nucleation of the new phase leads to a redistribution of mass in space as the phases 
have different densities in general. For the case D — 3 the function f(R) has been calculated in Ref. 0, where it 
has been assumed that both phases are incompressible quantum liquids. For the case D = 2 this approach gives 
a divergent kinetic energy and, consequently, we need a more detailed description of the phases 1 and 2 to reach a 
sensible result. Similarly, the effective mass per area of the boundary between the two phases cannot be calculated 
assuming incompressible phases. 



9 



Briefly summarizing, we have analyzed the problem of depinning of a (d + n)-dimensional massive elastic manifold 
from an 0(n)-invariant trapping potential V(u) = V(|u|) in the presence of a small external force F. For the case 
1 < d < 6 the Euclidean action has been calculated in the whole temperature range, see Eq. ([l4|) for T < T c (T c 
is given by Eq. (^)) and Eq. ( |l5| ) for T > T c , and we have found a sharp transition from quantum to classical 
behavior, see Fig. 3. The high-temperature asymptotics for the preexponential factor of the (1 + l)-problem has been 
calculated. Possible applications to the thermal depinning of vortices in high-T c superconductors and the nucleation 
problem in the vicinity of a boundary in D-dimcnsional phase transitions have been discussed and we have shown 
that for D = 2 — 4 nucleation at the boundary is more favorable than from the bulk. 
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